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Abstract —The handover process is one of the most critical 
fnnctions in a cellular network, and is in charge of maintaining 
seamless connectivity of nser equipments (UEs) across mnitiple 
cells. It is usnally based on signal measnrements from the 
neighboring base stations (BSs), and it is adversely affected by the 
time and frequency selectivity of the radio propagation channel. 
In this paper, we introdnce a new model for analyzing handover 
performance in heterogeneous networks (HetNets) as a function 
of vehlcnlar user velocity, cell size, and mobility management 
parameters. In order to investigate the Impact of shadowing and 
fading on handover performance, we extract relevant statistics 
obtained from a 3rd Generation Partnership Project (3GPP)- 
compliant HetNet simulator, and subsequently, we integrate these 
statistics into onr analytical model to analyze handover failure 
probability under flnctnating channel conditions. Computer sim¬ 
ulations validate the analytical findings, which show that fading 
can significantly degrade the handover performance in HetNets 
with vehicular users. 

Index Terms —3GPP, Bertrand’s paradox, handover, heteroge¬ 
neous network, long term evolution (LTE), mobility management, 
radio link fallnre, time to trigger. 

I. Introduction 

The new generation of wireless user equipments (UEs) have 
made user data traffic and network load to increase in an 
exponential manner, straining current cellular networks to a 
breaking point |[^. Heterogeneous networks (HetNets) which 
consist of traditional macrocells overlaid with small cells 
(e.g., picocells, femtocells, phantom cells etc.), have shown 
to be a promising solution to cope with this wireless capacity 
crunch problem 0. Due to their promising characteristics, 
HetNets have gained much momentum in the wireless industry 
and research community during the past several years. For 
instance, there have been dedicated study and work items 
in the third generation partnership project (3GPP) related to 
HetNet deployments 0 - Their evolutions are also one of the 
major technology components that are being considered for 
5G wireless systems Q. 

Despite their promising features, HetNets have introduced 
new challenges, such as the mobility management. Handover 
among different base stations (BSs) is the main process that 
supports seamless connectivity of UEs to the network. Due 
to increased number of cells in the network, it is difficult to 
support seamless mobility of UEs in a HetNet scenario since 


handovers may fail Q. In particular, using the same set of 
handover parameters of a traditional macrocellular network for 
a HetNet scenario will degrade the mobility performance of the 
UEs Q. For handover in Long Term Evolution (LTE) systems, 
signal measurements obtained at a UE from the neighboring 
BSs are reported by the UE to its serving BS, and the handover 
decision is made by the serving BS. In LTE HetNets, due to 
the small cell sizes, such measurement reporting by the UE 
may not be hnalized sufficiently quickly, and this might result 
in severe handover failure (HE) problems for the high velocity 
users Q. 

Handover performance has been studied for homoge¬ 
neous @-||T4| and heterogeneous i|Tg-@ network deploy¬ 
ments in the literature. In homogeneous networks, the authors 
in Q use computer simulations to investigate the handover 
performance of LTE networks, considering different measure¬ 
ment hltering parameters at the UE. A novel self-organizing 
handover management technique is proposed in ©-HD. 
where the network autonomously configures the mobility man¬ 
agement parameters for different scenarios, thereby improv¬ 
ing the handover performance of the homogeneous cellular 
network. Handover parameters (e.g. time-to-trigger (TTT), 
hysteresis threshold, etc.) are optimized in G3 to achieve 
robust and seamless mobility of UEs in a HetNet scenario. 
In 116) , mobility performance of UEs is evaluated in the 
co-channel small cell networks scenario; when the density 
of the small cell increases, switching off the macro cell 
is shown to provide seamless mobility for the low speed 
UEs, while it degrades the handover performance for the 
high speed UEs GD- Furthermore, in GD authors show that 
using intercell interference coordination (ICIC) techniques can 
enhance the handover performance for both low and high 
speed UEs. Mobility state estimation is performed in | fT9l to 
estimate the velocity of the UEs and thereby keeping the high 
speed UEs to macrocells and low speed UEs are offloaded 
to picocells, thereby enhancing the handover performance 
of the UEs. In | [20) , mobility performance is analyzed with 
and without inter-site carrier aggregation for macro and pico 
cells deployed on a different carrier frequencies. The authors 
in 0, @ aim to improve the mobility performance of UEs 
across different network types such as WiFi, WiMAX, LTE, 
and Bluetooth, by performing a vertical handoff (VHO). 
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Despite all these related work on mobility management in 
HetNets, there are only limited theoretical studies that ana¬ 
lyze the handover performance in HetNet scenarios. In 
the authors derive the handover rate and sojourn time of a 
UE for the Poisson-voronoi and hexagon cellular topologies. 
Expressions for call block and drop probabilities in a small 
cell scenario are derived in | |25) . Theoretical analysis for 
handover performance optimization is done in to quantify 
the user performance as a function of user mobility parameters. 
In a mathematical framework was proposed to model 

the handover measurement function, and expressions were 
derived for measurement failure and best target cell. In | [28) , 
handover performance analysis was performed as a function of 
handover parameters, e.g. TTT and UE velocity, and in |29| 
the analysis was extended to consider layer-3 measurement 
hltering process at the UE. To the best of our knowledge, apart 
from our preliminary results in pO) , there are no analytical 
results in the literature that study the HE probability in HetNets 
as a function of different mobility management parameters and 
under fading channel conditions. 

The main goal of this paper is to introduce a simple yet 
effective model for analyzing handover failures in small cell 
deployments, considering all important mobility management 
parameters of interest. Handover trigger locations at a picocell, 
and radio link failure locations at a macrocell and a picocell 
are modeled using co-centric circles. Considering a linear 
mobility model for UEs, HE probabilities for macrocell and 
picocell UEs are derived in closed form for various scenarios. 
The analysis is then extended to fast-fading and shadowing 
scenarios; relevant statistics in a fading scenario are extracted 
from a 3GPP compliant system level simulator, to facilitate 
semi-analytic expressions for handover failure probabilities. 
All theoretical results are validated through simulations, where 
impact of different parameters on handover failure and ping- 
pong probabilities are investigated. 

The paper is organized as follows. In Section the han¬ 
dover process in LTE and the handover measurement process 
in a UE are reviewed. In Section a new analytical model 


for handover performance analysis is presented. In Section IV 


HE probability expressions are derived in the absence of fading 
and shadowing, while a semi analytical approach for analyzing 
the handover performance in fast fading environments is intro¬ 
duced in Section |V] In Section |Vl] the theoretical expressions 
for the HE are verihed via simulations, and the last section 
provides the concluding remarks. 


H. Review of the Handover Process in LTE 
A. Different Stages of Handover Process 

The key steps of a typical handover process in a HetNet 
scenario are illustrated in Eig. [T] Handover decisions are 
based on the signal strength measurements of the neighboring 
BSs done at the UE. In LTE, UEs perform reference signal 
received power (RSRP) measurements to assess the proximity 
of neighbouring cells. An example for the downlink (DL) 
RSRP measurement prohle of a macrocell and a picocell, 
measured by a mobile UE, are shown in Eig. [T] Once the 
measurements are performed, the UE checks for the handover 


event entry condition, e.g., when the signal strength Pp from 
target cell (e.g., a picocell) is larger than the signal strength 
from a serving cell (e.g., a macrocell) Pui plus a hysteresis 
threshold (step-1). Even when this condition is satished, the 
UE waits for a duration of TTT, before sending a measurement 
report to its serving cell (step-2). 

The use of a TTT is critical to ensure that ping-pongs 
(successive and unnecessary handovers among neighboring 
cells), generated due to fluctuations in the link qualities from 
different cells, are minimized. If the handover event entry 
condition is still satished after the TTT, the UE sends the 
measurement report to its serving BS in its uplink (step-3), 
which then communicates with the target cell. If both cells 
have an agreement and the handover is to be performed, the 
serving BS sends a handover command to the UE in the DL to 
indicate that it is should connect to the target cell (step-4). The 
handover process is hnalized when the UE sends a handover 
complete to the target cell, indicating that the handover process 
was completed successfully (step-5). 



Fig. 1: Handover failure problem in HetNets due to small cell size. 


B. Handover Measurements Procedure in LTE. 


Different measurements obtained at a UE during a han¬ 
dover process are summarized in Eig. It is important to 
note that the RSRP measurements and Pp at a UE are 
obtained after a hltering process for smoothing the signals, 
in order to mitigate the effects of channel huctuations. The 
hltering of handover measurement is performed in Layer- 
1 (LI) and Layer-3 (L3), as shown in Eig. m Initially, 
the UE obtains each RSRP sample by linear averaging over 
the power contribution of all reference symbols carrying the 
common reference signal within one subframe (i.e., 1ms) 
and also over the measurement bandwidth of at least six 
physical resource blocks. In a typical handover measurement 
conhguration shown in Pig. |^[b) the UE performs the layer- 
1 hltering by obtaining an RSRP sample every 40 ms, and 
performs linear averaging over a number of successive RSRP 
samples, usually 5 samples. As a result, LI hltering performs 
averaging over every 200 ms to obtain an LI sample, Min), 
given by |j^ 


1 


M (n) = - RSRPLi{5n — k) , 


( 1 ) 


fc =0 
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(b) Processing of the RSRP measurements through LI and L3 
filtering at a UE 0. 

Fig. 2: Handover measurement performed by the UE through filtering 
process in LTE. 



where n is the discrete time index of the RSRP sample, 
RSRPi^i is the RSRP sample measured at every 40 ms by 
the UE, and k is the delay index of the filter. 

A UE further averages the LI samples through first-order 
infinite impulse response filter called L3 filtering, which is 
given by 

F{n) = (1 - a)F{n - 1) -I- al01ogio[-W(n)] , (2) 

where a is the L3 filter coefficient. 

Then, UE periodically checks the resulting L3 sample every 
Td seconds (e.g., 200 ms in 3GPP LTE | [T^ ) for the handover 
entry conditions. Td refers to L3 sampling period in the 3GPP 
LTE standard. If the handover entry condition is satisfied, then 
rest of the handover steps may follow as described previously. 

III. A Geometric Model for the Analysis of 
Handover Failures 

In order to evaluate the handover performance in HetNets, 
there is a standard hotspot model with specific simulation 
scenarios and parameters presented in the 3GPP study item 0. 
This hotspot model is based on a bouncing (hotspot) circle 
concentric within the picocell, whose radius is assumed to be 
200 m. The starting UE position is chosen randomly on the 
bouncing circle and the UE follows a linear trajectory. The 
UE does not change the direction until it hits the bouncing 
circle and bounces back with a random angle, as shown in 
Fig . 1^ With such a model, theoretical analysis of the handover 
failure probabilities is challenging, due to the complexity of 
modeling the statistics of a UE’s sojourn time with in a 
picocell. Instead, we propose a simple geometric model in the 
following subsection for the handover performance analysis. 


A. New Geometric Model for Handover Performance Analysis 

In order to simplify the handover model shown in Fig. 
we initially consider the handover metrics of a single user 
in the absence of fading, and develop a new framework to 
facilitate closed form analysis of HE probabilities. Later, using 
this framework as a reference, we extend our analysis of HE 
probabilities into the scenario where there is channel fading. 

Considering the handover measurement process summarized 
in Fig. 1^ and assuming that there is no fading, we investigate 
where the TTT are initiated within the coverage area of a 
picocell and also the handover failure locations experience 
by the UE. Using a 3GPP compliant simulator and based 
on the assumptions shown in Fig. 0 0’ handover trigger 
(red dots, after step-1 in Fig. and HE locations (blue 
cross) are obtained. Fig. shows the ring of handover trigger 
locations as a result of discrete measurement carried out at 
UE. These handover trigger locations extend inwards from the 
ideal coverage area of a pico BS (PBS), since a UE may delay 
initiation of the TTT due to the filtered RSRP measurements 
being available only with 200 ms intervals. Fig. also shows 
that if we neglect the impact of sectorized cell structure, the 
HE locations (the blue cross signs where the wideband signal 
to interference plus noise ratio (SINR) becomes lower than a 
threshold) can be approximated by a circle. 

Based on Fig. we model the picocell coverage and HE 
locations geometrically as concentric circles shown in Fig. 
with radius R and respectively. The HE locations shown 
in Fig. are for macro to pico HFs and in the same way 
it is reasonable to approximate pico to macro HE locations 
as another concentric circle with r-p shown in Fig. In 
Fig. 0 V denotes UE velocity, 9 is the angle of UE trajectory 
with respect to horizontal axis, T^, Tp are the TTT duration 
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for MUE and PUE respectively. In the following subsection, 
we incorporate the discrete measurements process in our 
geometric model by proposing a standard distribution that fits 
to the handover trigger locations. In the subsequent subsection, 
we use it to model the UE’s sojourn time in order to facilitate 
the theoretical analysis of macro-cell UE (MUE) and pico-cell 
UE (PUE) HE. 
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Fig. 4: Flandover trigger and handover failure locations in an example 
picocell using a 3GPP-compliant simulator. The MBS is lo cate d at 
(1500,1500) m and has three sectors (Tj = 200 ms) |28|. 



Fig. 5: Geometric handover model to analyze HFs for MUEs and PUEs, and 
considers the effects of L3 filtering. 


B. Modeling the Handover Trigger Locations 

To model the statistics of the handover trigger locations 
which are offset due to L1/L3 filtering, we examine the 
distance of each handover trigger location shown in Eig. 
from the ideal picocell coverage boundary, which we define 
as handover offset distance. In this paper, we model the 
distribution of the handover offset distance using a uniform 
distribution for the scenario with no fading, and using a custom 


distribution for the scenario with fading (to be studied in 
Section [V|. 

We incorporate the discrete measurement process in our 
geometric model through modeling the handover trigger loca¬ 
tions which is explained Eig. Let us consider a UE traveling 
with a velocity v. It is assumed that the UE is traveling in a 
straight line and due to the discrete measurements performed 
by the UE every seconds, it checks whether the handover 
condition is satisfied at integer multiples of T^. Therefore, 
when the UE crosses the picocell coverage circle, TTT is 
not initiated immediately. It has to wait for the remaining 
duration of time seconds depending on where the UE 
processed the measurement before crossing into the picocell 
coverage area. Let us consider different starting instance of 
UE crossing the picocell coverage area shown in Eig. 0 In 
the first instance, the UE starts as an MUE at green dot and 
travels through the picocell coverage circle with a distance 
vT^. We can see that the TTT is triggered after the distance 
vT^. In the second instance the MUE starts earlier than the first 
instance. Therefore, after a distance vT^, the location where 
the TTT would be triggered after L1/L3 filtering will be closer 
to the ideal picocell coverage boundary than the first instance. 
In this paper, we model the distance between the triggering 
location of the TTT and the picocell coverage circle as a 
random variable, and we denote this distance as rd. If we 
consider all possible instances, we assume that the distance 
from the cell edge reference point is uniformly distributed, i.e. 
rd ~ if [0, tiTd), which models the handover trigger locations. 

To verify our model, we check the histograms for handover 
offset distance rd for the trigger locations generated by the 
3GPP-compliant simulator. Based on the assumptions shown 
in Pig. and handover parameters shown in Table |Ij the 
handover trigger locations are generated for the following three 
cases: 

1) Without shadowing and fast-fading (Case-1), 

2) With shadowing, but without fast-fading (Case-2), 

3) With shadowing and fast-fading (Case-3). 


TABLE I: Handover pai'ameter sets. 


Profile 

Set-1 

Set-2 

Set-3 

Set-4 

TTT (ms) 

480 

160 

80 

40 

L3 filter coefficient (a) 

4 

1 

1 

0 


The handover offset distance histograms for Case-1 is shown 
in Pig. 1^ and we can see that for higher UE velocity the 
probability density function (PDP) flattens according to our 
proposed uniform discrete assumption, as described earlier. 
Modeling the handover process in the fading scenario follows 
from Case-1, and is achieved by finding the handover offset 
distance histograms for the trigger locations in the Case-2 and 
Case-3, which is studied in Section |V] 


C. Modeling the UEs’ Sojourn Times 

The sojourn time estimated in a picocell may be different 
when using the geometric model in Fig. or when using the 
bouncing ring model of Fig. In order to justify the use of our 
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Fig. 6: Handover offset distance histograms for i; = 60 km/hr in Case-1. 
For the scenario of = 200 ms adopted in 3GPP LTE, the handover offset 
distance histogram can be reasonably modeled using a uniform distribution. 


geometric model in Fig. to model the scenario in Fig. we 
consider the Bertrand’s Paradox and the three probability 
density functions described therein. We will show that one of 
these probability density functions well matches the behaviour 
of the bouncing ring model, thus validating our modelling. 

The Bertrand’s Paradox studies the probability that a ran¬ 
dom chord of a circle with a radius R is larger than a threshold 
parameter. In essence, this probability leads us to the statistics 
of the sojourn time in a given picocell. Due to different 
interpretation of randomness of a chord in a circle, there are 
three different models for the PDF of the chord length. 

Model 1 : When we choose randomly two points on 
a circle and draw the chord joining them, without loss of 
generality, we may position ourselves at one of them and 
examine the relative location of the other points. If angle 6 
is uniformly distributed between [—f, f ], then the PDF of the 
chord length I is given by: 


/i(0 


2 

7rV4i?2 _ 12 ■ 


(3) 


Note that this interpretation corresponds to the model used in 
Fig. ID 

Model 2 : If we choose a chord whose direction is fixed 
and perpendicular to a given diameter of the circle, then we 
assume that the point of intersection of the chord with the 
diameter has a uniform distribution. Therefore we can assume 
that the perpendicular distance r = B? — ^ from the chord 
to the center of the circle is uniformly distributed between 
The PDF of the chord length I is then given by: 


/ 2(0 


I 

2i?v/4i?2 - l^ ■ 


(4) 


Model 3 : A chord is uniquely determined by its mid¬ 
point, for which a perpendicular line extending from the circle 
center intersects with the chord. We assume this intersection 
point is uniformly distributed over the entire circle and PDF 


Fig. 7: PDFs of the chord length for the three different approaches in 
Bertrand’s Paradox, and the histogram of the chord length for the 
bouncing-ring simulations in Fig. 


of chord length is given by: 


Ml) 


i 

^ ■ 


(5) 


In order to evaluate how closely the three approaches in the 
Bertrand’s Paradox capture the picocell sojourn time in Fig. 
we compare the PDFs of chord lengths for the three Bertrand’s 
Paradox cases, with the simulated chord length histogram for 
the bouncing ring model in Fig. Considering R = 21.7m 
and plotting the histograms of chord length overlayed with 
PDFs of all the three solutions in ([^-Q, we obtain the 
results in Fig. Model 1 shows a reasonable match with 
the histograms and therefore we use the PDF given in ([^ for 
our MUE and PUE HE analysis in both no fading and fading 
scenarios. 


IV. Handover Eailure Analysis without Eading 

Using the geometric model in Eig. |^and the PDE of chord 
length in ([^, in this section, we derive the handover failure 
probabilities for MUEs and PUEs considering L3 filtering 
for the ideal handover case. We consider that a UE checks 
the handover entry condition at every sampling period of 
the L3 hlter. When the handover condition is satisfied, then 
TTT of duration is triggered. The rd is a random variable 
accounting for discrete measurement interval carried out in the 
UE through L3 filtering, and we assume that rd is uniformly 
distributed between [0,t;Td], yielding the following PDE 


fM) = 



0 < Td < uTd 
otherwise 


(6) 


In other words, rd models the random offset between 
the intersection of the UE trajectory with the ideal picocell 
coverage circle, and the location when the filtered measure¬ 
ments become available to the UE after entering the picocelTs 
coverage area. In the following section, we derive the no 
handover (NHO) probability for MUEs and HE probabilities 
for both MUEs and PUEs using Eig. and equations 0 
and (|^. 
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A. No Handover Probability for MUEs 

After TTT of duration Tm is triggered, the MUE does not 
make a handover if it leaves the picocell coverage circle before 
the end of T^. In Fig. |^we can see that the MUE does not 
make a handover if vT^ + rd is larger than the chord length 
l{9) = 2i?cos(0), and for which the chord does not intersect 
with the MUE HE circle. Depending on the value of vT^ 
relative to rd and picocell coverage, NHO probability for a 
given UE can be analyzed for three different cases. 

1) vT^ < 2^R? — r^ — vT^\ Since random variable rd 
and 9 are independent, we can just multiply the two PDFs 
in 0 and 0 to obtain the joint PDF. Then, after some 
manipulation, NHO probability is expressed as 


3) If UE velocity is too high and product vT^ is itself 
very large and greater than the chord length 2-y/i?^ — r^, then 
there is no effect of L3 sampling period on NHO probability. 
In this case, we obtain a constant NHO probability due to the 
hxed limits in 9, and using 0 the NHO probability becomes 


2y/R^-rl 

P„„o = P(i(») > 2,/W^) = I 

0 



( 11 ) 


*^NHO — 


< vT^ + rd) 




vTd i>r„,+rd 

f 1 f 


f dr drd , 


7rV4i?2 - P J vTj 7rV4i?2 _ p 

0 uTm 

V4i?2 - {vTdd, + vTaY + {vT^ + t;Td)fd 


TTuTd 

- ^/AR^ - {vT^Y - vT^ tan-i 


vT„ 


V4i?2 _ 


where 


Td = tan 


-1 


uTm + vTd 


y4i?2 - [vT^ + vTdf 


(7) 


( 8 ) 


2) vTtn > 2^R? - r^ - vT^. If vT^n + vT^ is greater 
than the chord length 2y/RF^^r^, then we have to subtract 
the integration region which overlaps with the MUE HE circle 
from case 1, which yields 


B. HE Probability for MUEs 

When MUE reaches the MUE HE circle before TTT expires, 
there will be handover failure and the MUE fails to send 
a measurement report. Therefore, the MUE fails to connect 
with the pico cell. From Fig. we can see that if l{9) > 
2dJR'^ — r^, then the MUE trajectory intersects with the MUE 
HE circle and therefore there may be a MUE handover failure. 

A MUE HE occurs only when the total distance traveled 
by MUE, which is vT^n + ^d, is greater than the distance 
dHF,m(^, 7?, Tin). This distance refers to total distance travelled 
by the UE from the ideal picocell coverage to the MUE HE 
circle and is given by 


dHF,m(6', 7?,rm) = Rcos{9) - \Jr'^- R?s\t?{9) . (12) 


To obtain the MUE HE probability, hrst we evaluate the MUE 
HE condition vT^ + rd > dHF,m(7*, 7?, rm) in terms of UE 


trajectory l{9). Using (12i, we can write 


T^nho = P(z(0) < wTjn + rd I rd < 2 a/ R? - - vT^ 

vTm uTd 'uT'm + T'd 

= /—7=^=dZ+ f - ^ - drdrd 

J 7r/4i?2 -p J vTd J 7rV47?2 - p 

0 0 vTm 

vTd uTm+uTd 

- — drdrd ■ (9) 

J vTd J 7rV47?2 - p 

The limits of the random variable rd in the third term represent 
the integration region which overlaps with MUE HE circle. 
After solving the integral in 0 and using Td we get 


vTjn + rd > 7?cos(a) — \Jr'^ — R^ sin^(a) . (13) 

Using l{9) = 2Rcos{9), we get 9 = cos“^(^^), and using 
this in we get the MUE HE condition in terms of l{6) 
as 

m> +(vT^ + ra) . (14) 

w7m + rd 

Then the MUE HE probability is calculated differently for the 
following four cases. 


T^nho = 


TTwTd 
— vTjn tan“^ 


ivT„, + vTd + 2s/R^ - r^)Td - 2r„ 
vTm 

V4i?2 - (vT^) 


+2y/AR^-{vT,n + vTd)^ 




(2v^T?2 - r^)tan ^ 



1) vTdn < R — Tm — vTd'. We know that for 0 = 0, 
dHF.min = R — fm- If the Velocity of the UE is very slow then 
the distance traveled by UE including the discrete L3 hlter 
measurement offset is vT^a + rd. If this distance is less than 
dHF.min, then there will be no MUE HE, i.e., T^HF.m = 0. 

2) R- < vT^ + uTd < s/W^TfT. when the 
velocity of the UE is high, then we have obtained the MUE 
HE condition in ([T4| and we obtain the MUE HE probability 
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as 


Pt 


HF,m — 


»(Z(0) 


> 


' m 

+ Td 


iiTd 




+ + rd)) 

9 

dZ drd 


0 

2R 


+ 


+t;T„ +rd 7r\/4i?2 _ ;2 


d/ drd 


7rV4i?2 - P 


r'^^d 


= 1 - 


irvTd 


/i(rd)drd , 


where /i(rd) is 


tan 


-1 




(15) 


v/4i?2(r;r^ + rd)2 - {R^ - + (vT^ + , 

(16) 


3) vT^ > — r'^ — vTd'. For this case, the MUE 

HE probability is the same as for case 2. 

4) If the UE velocity is too high and is very large 

and greater than the chord length then MUE HE 

probability is independent of the random variable rd, and we 
obtain a constant probability, which is the same as for case 3 
of NHO probability. 


C. HF Probability for PUEs 

In order to observe a PUE HE there should be a successful 
handover of MUE to the pico-cell. After a successful handover, 
the PUE starts recording the measurements of PNB for every 
L3 sampling period T^. If a PUE enters the coverage of the 
macrocell and handover event entry condition is satsified (e.g., 
L3 filtered RSRP of the macro cell is larger than that of the 
picocell plus an hysteresis), then TXT of duration Tp (see 
Fig. is triggered. For simplicity, we assume that discrete 
offset random variable rd is the same whenever there is a 
handover to picocell or macrocell. 

If a PUE reaches the PUE HE circle before the TXT expires, 
there will be a PUE HE. In other words, a PUE HE occurs 
when the total distance travelled by PUE {vTp + rd) is greater 
than the distance d}iF^p{9, R,rp). If we consider a point on 
the ideal picocell coverage area where the UE starts entering 
the coverage of the macrocell after a successful handover to 
picocell, then distance from this point to the PUE HE circle 
is given by 

rfHF,p(0, R, fp) = Rcos{6) + \Jrp — R? sin^(0) — d{9) . 

(17) 


Xo obtain the PUE HE probability, we will evaluate the 
PUE HE condition vTp + rd > dYi-p,p{9,R,rp) in terms of 
UE trajectory l{9) like we did before for MUE HE. Using 
0’ we get the condition for observing PUE HE as 

m>^y^-{vTp + rd) . (18) 

vTp + rd 


Based on the condition in (18i and the condition l{9) > 
R^ — r^, we can show that when uTm + vTp > 


\J'^p — r^ — yR? — r^, there will be a PUE HE. Eor different 
values of vT^ and vT^, the PUE HE is given as follows: 

1) vT^ + vTd < y/R"^ — r‘^{- In this case, in order to 
observe PUE HE, we have to make sure that there is no MUE 
handover failure and we will use MUE HE condition in ( [l4| ) 
and PUE HE condition in ( fTS) ) to obtain PUE HE. We can 
see that there is no MUE HE even when l{9) > 21^7^2^77^. 
Moreover, there will be a possible PUE HE for the condition 

vT^ + vTp > ^rl - r^ - ^- r^ . (19) 


After some manipulation PUE HE probability is given by 

rl - i?2 


Fhf.p — ■ 


max ( vT^ + rd, - 

vlp + rd 


- {vTp + rd)^ 


< d{a) < min 


. / i?2 - ' 


\vT„, 


rd 


+ ('L’T'ni + rd), 2i? 


Consider the following definitions for brevity: 
r^ — R^ 

dvp = - {vTp + rd) , 

vTp + rd 

j^2 — j,2 

drm = ^ " + {vT^ + rd) 

vTra + rd 

^vm — 4” rd . 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


Using (21 h23i we can then write Lp = max(dvm, ^rp) and 
Lin = niin(drni, 2R). Using the PDE in ([^, we can calculate 
the PUE HE probability after some manipulation as 




Jo nViR^ - R 


R—r^ — vTj^ 


vTd 


0 

— J tan” 
0 


nvTd 

R—Tra — vTm 


4i?2 - 42^ 


tan 


d^ drd 

drm 


V4i?2 - < 


drd 


vTd 

drd — / tan 


-1 ( 


\yAR^-di 

(24) 


where lp is 

y_8i?2 + 8r2 + (t;Tin)2 - 2v^T^Tp + (vTp)^ - ^v^TddTp 

4 

(25) 


2) yR-^ - r^ < vT^ + vTd < 2i/fl2 - r^: In this 
case, MUE HE occurs when l{9) > 2^/772 — Using this 
condition and the PUE HE condition in ( [T8] l, the PUE HE 
probability is given by 

/ / 7*^ — \ 

Fhf,p = P ( max + rd, - (r^Tp + rd) j 

< dia) < 2yR? - r2^^ . 

( 26 ) 


drd , 
















































Using Q, ([2T]|-(23 I, Lm and Lp, we can find PUE HE 
probability after some manipulation as 

rvT^ 1 2 


= — tan 

TT 


-1 


vTa 


tan 


-1 


v/4i?2 - dl. 


+ J tan 
0 


drd . 


-1 




drd 


(27) 


5) t;r„i + vTd > 2^i?2 - r^: In this case, PUE HE is 
the same as in case 2 of PUE HE probability. 

V. Handover Failure Analysis with Fading 


The geometrical model framework to perform handover 
failure analysis in the fading scenario is similar compared to 
no fading scenario. We model the handover failure locations 
in the fading scenario as a circle, because the wideband SINR 
(which dictates handover failures) is typically averaged over 
a large bandwidth, and the effect of fading is minimal. On 
the contrary, the handover trigger locations in the fading 
scenario depend on the RSRP, which is measured over a 
narrow bandwidth of six resource blocks in LTE. Therefore, 
RSRP is subject to much larger randomness, and hence, 
may significantly affect the handover trigger locations. In 
the following section, we explain the modeling of handover 
trigger location in the fading scenario and investigate its 
distribution to facilitate the theoretical analysis of handover 
failure probability. 


A. Modeling the Handover Trigger Locations in the Fading 
Scenario 

Modeling the handover trigger locations in the fading sce¬ 
nario is performed by examining the handover trigger locations 
shown in Fig. Comparing Fig. with Fig. we can 

1400p -,-,- 111,1 ^ 

• Handover trigger locations in fading scenario 
O Picocell coverage circle center 


1350F 



11001 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

800 850 900 950 1000 1050 1100 1150 1200 

X-axis (meters) 

Fig. 8: Handover locations in the fast-fading and shadowing scenario 
(Case 3) for UE velocity i; = 60 km/h and = 40 ms. 

see that due to fading channel conditions handover trigger 
locations in Fig. are not delimited within the coverage 


boundaries of the picocell. In other words, TTT timer can be 
initiated for locations that are far away from the ideal picocell 
coverage area. To model the statistics of the handover trigger 
locations with L1/L3 filtering, we examine the distance of 
each handover trigger location shown in Fig. from the ideal 
picocell coverage boundary which we define as handover offset 
distance. We model this distance using a random variable 
denoted by and obtain histograms for it as shown in 
Fig. 0 The negative distance in histograms is due to the 
possibility of handover trigger locations being outside the ideal 



Distance (meters) 

Fig. 9: Handover offset distance histograms for v = 120 km/h in fading 
scenario. 


standard distributions, we conclude that there is no standard 
distribution that reasonably approximates the histograms in 
Fig. 0 Therefore, for the fading scenario, we directly use the 
histogram data to obtain the semi-analytic handover failure 
expressions for MUEs and PUEs. The histograms for fd is 
given by /(fd), where fd € (f^in, Tmax)- 


B. MUE HF Probability Analysis for MUEs 

The derivation of MUE HF probabilities in shadowing 
and fast-fading is carried in a similar manner to that of the 
ideal handover model case. The only difference comes with 
distribution of the random variable fd, which is obtained 


directly from the histograms as discussed in Section V-A The 


MUE HF probability for a given UE can be analyzed as 

2 2 ^max 

^ + ivT^ + fd)) = f /(fd)dfd 

V vT^ + rd / J 


2R 


fmin fmin 

' max V m 

= J f{fd)dfd- J Ii{fd)f{fd)dfd 


R—rjy,—vTjj. 


R—Tra—vTm 


-J tan ^ ^ -^/(Td)dTd - J f{fd)drd , (28) 

R'^-rf-vT^ 





































9 


where /i(rd) is 


tan 


-1 


-rl^ + {vTra + fdf 


+ fd)2 - (i?2 _ + f^)2)2^ 

(29) 


In Section VI the above integrals are solved numerically to 


provide numerical results for handover failure in the presence 
of fading. 


C. PUE HE Probability Analysis for PUEs 

The derivation of PUE HE probabilities in shadowing and 
fast-fading is carried out in a similar manner to that of the 
ideal handover model. The only difference is in the distribution 
of the random variable fd, which is obtained directly from 


the histograms as discussed in Section V-A Using MUE HE 
condition in ( [l4| ), PUE HE condition in (|18|), and after some 
manipulation, the PUE HE probability is expressed as 


-Phf.p = ■ 


r^ — R? 

max ( uTin -f fd, ^- (vT-p + fd) 


vTp + Td 


< d(a) < min 


/ i?2-- 


\vT„, 




+ {vTra + fd),2R 


Consider the following definitions for brevity: 
— R? 

“ TTE ^d) , 


vTp + fd 


_ 2 

drm = ^ + {vT^ + rd) 

vT^ + rd 

^vm — -f fd ■ 


(30) 

(31) 

(32) 

(33) 


Then, using (31 h33i we can write Lp = max((ivm, (^rp) 

Lm = min((iri„, 2R). Using the PDE in we can calculate 
the PUE HE probability as 


^ max -^m 


R—rm — vTm 


r^ — R—vT^ 


J /(^d)dfd-y /(fd)dfd 


2R-vT^ 

— / tan" 


2R-vTm 


Ip 

— J tan 

r^ — R—vTr 


-1 


V4i?2 - dl 

dr 




4i?2 - d^p 


dfd 


dfd 


+/‘“" 


/(rd)dfd 


rjR^-rf-vTrr, 

r 2 


J TTvTd 

R—rm — vTm 


tan 


-1 


4ii2 _ ^2 


dfd . (34) 


Again, to obtain numerical results in Section |VI| the above 
integrals are solved numerically using the histograms /(r^) 
from Eig. 


VI. Numerical Results 

In order to validate our analysis, computer simulations are 
carried out in Matlab. Initially, as shown in Eig we fix the 
starting position of MUE to a reference point on the picocell 
coverage circle. MUE travels a distance equal to uT^+rd from 
the reference point with different realizations of angle 6, which 
is uniformly distributed in Since rd is uniformly 

distributed between [OjUTd], we have different values of rd 
for each 0. After MUE travels a distance equal to vT^ + rd, 
the final point of the MUE is checked for intersection with 
MUE HE circle. If it intersects with the MUE HE circle then 
there is a MUE HE. We aggregate all these MUE HE events 
and normalize over 0 and rd to obtain MUE HE probabilities 
for each UE velocity. 

In order to obtain PUE HE probability, we find end points 
of the chord using l{9) — 2iicos(0) for the picocell coverage 
circle for different 9 realizations. Then, we take those points as 
the reference point and fix the starting position for the PUE. 
Next the PUE travels a distance equal to vTp + rd and the 
final point of the PUE is checked for intersection with PUE 
HE circle. If it intersects with the PUE HE circle then there is 
a PUE HE. We aggregate all these PUE HE events and divide 
with the number of 9 and rd realizations to calculate PUE HE 
probabilities. In the following subsections the MUE and PUE 
HE probabilities are shown for ii = 64 m, rm = 50 m and 
rp = 78 m in the no fading and fading scenarios. 


A. Results with No Fading 

Using the above simulation assumptions, theoretical MUE 


HE and PUE HE probabilities derived in Section IV 


plotted as a function of UE velocity and are verified via 
simulation results. Initially, it is assumed there is no fading 
or shadowing. Then, the MUE HE probability is shown in 
Eig. lQ[a) for different TTT and Td values. We see that as UE 
velocity increases, the MUE HE probability increases. On the 
other hand, when the sampling period of L3 filter decreases, 
the MUE HE probability decreases, since the TTT can be 
initiated earlier. Eor example, the MUE HE probability for UE 
velocity 80 km/hr improves by approximately 10 percent when 
sampling period is reduced from 150 ms to 50 ms. By reducing 
TTT to 160 ms we can see that the MUE HE probability 
becomes almost negligible. This is because the UEs finalize 
the handover in a quicker way, when compared to larger TTT 
values. 

PUE HE for different TTT and Td values are shown in 
Eig. iq[b) We can see that as UE velocity increases, PUE 
HE probability increases and it is improved when sampling 
period of LI and L3 filter decreases. Eor example, the PUE 
HE probability for UE velocity 120 km/hr improves by ap¬ 
proximately 5 percent when sampling period is reduced from 
150 ms to 50 ms. By reducing the TTT to 160 ms we can see 
that there will be no PUE HE probability when no fast fading 
and shadowing are considered in analysis/simulations. 
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(b) PUE HF probability with no fading. 

Fig. 10: Theoretical (lines) and simulated (markers) no fading results as a 
function of UE velocity for /? = 64 m, Vm = 50 m and rp = 78 m. 


The downside of reducing the sampling period of L1/L3 
filter is that it increases the unnecessary handovers called ping- 


pongs shown in Fig. 11 We can see that ping-pong probability 


increases when the sampling period of L1/L3 filter is reduced 
from 200 ms to 50 ms. The reason for this is that for a larger 
sampling period of L1/L3 filter, a UE will stay with its serving 
cell for a longer time before initiating the handover process, 
which will naturally reduce ping-pong handovers. 


B. Results with Fading 

We use the histogram data shown in Fig. to generate 
samples for fd in a fast-fading and shadowing scenario, follow 
the other simulation assumptions stated previously, and plot 
the MUE HF and PUE HF probabilities as a function of UE 
velocity. MUE HF probability plots are shown in Fig. 


We can see that MUE handover performance is degraded for 
all UE velocities compared to no-fading scenario in Fig. lC||^a) 


and it is improved when the sampling period (Td) of the filter is 
decreased. We can see that for MUE traveling with a velocity 
120 km/h, the MUE HF is improved by 12.19 percent when 
sampling period of the filter is reduced from 200 ms to 50 ms. 
Results show that even at low UE velocities, there may be 
on the order of 10% HF probability. Note that these results 
consider a worst-case simulation scenario, in which the UE 
starts its path at the coverage area of a picocell base station 
as shown in Fig. 

On the other hand, the PUE HF plots are shown in 
Fig. m[b)| We notice that the PUE HF probability increases as 
UE velocity increases and it is improved when the sampling 
period Td of the filter is decreased. The PUE HF probability 
improves by 3 percent for PUE traveling with a velocity 120 
km/h in the fading channel conditions. This is due to higher 
HF probabilities of MUEs at higher speeds. 

In order to investigate the impact of TTT and Td on ping¬ 
ping handover performance, ping-pong probability plots in 
the fading scenario are shown in Fig. 13 We can see that 
there are more ping-pongs when Td and TTT are reduced 
and this is because, using larger TTT and Td, UE will take 
more time to complete the handover process, as a result ping- 
pongs is reduced. In the Fig. [T^ we can see that there are 
more ping-pongs in fading scenario compared to no-fading 
scenario shown in Fig. 11 This is because, in the case of 
fading scenario and assuming from the picocell perspective, 
the the link quality of the picocell might decrease below the 
macrocell because of the channel imperfections. As a result 
there might be handovers back and forth and UEs time-of- 
stay will be less than the ping-pong threshold causing more 
ping-pongs. 


VH. Conclusion 

In this paper, we developed a simple geometric abstraction 
to derive analytic and semi-analytic expressions for MUE and 
PUE HF probabilities for scenarios with no-fading and fading, 
respectively. Despite the several simplifying assumptions con¬ 
sidered in the geometric model and framework, our findings 
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(a) MUE HE probability with fast-fading and shadowing (Case 3). 




(b) PUE HE probability with fast-fading and shadowing (Case 3). 


Fig. 12: Theoretical (lines) and simulated (markers) no fading results as a 
function of UE velocity for = 64 m, Vm = 50 m and rp = 78 m. 


still provide several useful insights. In particular, how different 
parameters explicitly affect the handover failure of UEs can 
be captured explicitly, which is not possible in the case of 
earlier simulation studies available in the literature. The results 
show that fading degrades the handover performance for all 
UE velocities. The handover performance for both MUEs 
and PUEs are improved when filtering sampling period is 
reduced. This improvement of handover failure performance is 
larger in the fading scenario. In particular, the handover failure 
performance improved by 12.16 percent for MUEs in the 
fading scenario, compared with the 3.92 percent improvement 
in the no fading scenario. The improvement in the handover 
failure performance due to shorter L3 sampling period comes 
at the cost of increased number of ping pong handovers, which 
will be studied further in our future work. 
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